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Excess of positrons in cosmic rays: A Lindbladian model of quantum electrodynamics 
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The fraction of positrons and electrons in cosmic rays recently observed on the International Space 
Station unveiled an unexpected excess of the positrons, undermining the current foundations of cos¬ 
mic rays sources. We provide a quantum electrodynamics phenomenological model explaining the 
observed data. This model incorporates electroproduction, in which cosmic ray electrons decelerat¬ 
ing in the interstellar medium emit photons that turn into electron-positron pairs. These findings 
not only advance our knowledge of cosmic ray physics, but also pave the way for computationally 
efficient formulations of quantum electrodynamics, critically needed in physics and chemistry. 
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Introduction. Cosmic ray electrons are primarily be¬ 
lieved to be generated form quasars, micro quasars, and 
supernovas Q]. Collisions of cosmic ray hadrons with the 
interstellar medium yield positrons, which are character¬ 
ized by a monotonic decrease with energy EHS|. How¬ 
ever, this explanation is challenged by recent measure¬ 
ments of cosmic ray electron and positron fluxes with 
energies ranging from 0.5 GeV to 1 TeV, detected by 
the Alpha Magnetic Spectrometer (AMS) on the Inter¬ 
national Space Station [TUTU] (see Fig. [l]). These ob¬ 
servations revealed an inexplicable rise of the positron 
fraction at higher energies, which is a subject of heated 
debate [MHma. In particular, an annihilation of dark 
matter has been entertained as the unaccounted source 
of positrons mmmm- Less speculative sources 
have also been put forth mnn], but these explanations 



FIG. 1: (color online) Comparison between the proposed 
model for electroproduction and the measured positron frac¬ 
tion in cosmic rays. 
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were found to be inconsistent with the observed boron- 
to-carbon ratio [ana- 

in this Article, we show that the positron excess can 
arise as a result of the cosmic ray electrons directly in¬ 
teracting with the interstellar medium. In fact, ultra- 
relativistic electrons injected into a medium produce 
electron-positron pairs through the process of electropro¬ 
duction mm- One possible channel of such a reac¬ 
tion is as follows: Under the influence of the external 
Coulomb field of a nucleus, the incident electron emits 
an energetic gamma quant that subsequently decays into 
a positron-electron pair. We propose a relativistic Lind¬ 
bladian model for electroproduction that not only quan¬ 
titatively describes the observed positron excess (see Fig. 
[I]), but also predicts a plateau for the positron fraction 
in the high energy limit. 

The Lindbladian model. In order to derive the mas¬ 
ter equation, we employ the paradigm of Operational 
Dynamic Modeling (ODM) [20f,' 2T], designed to deduce 
equations of motion from the dynamics of observable 
averages, which are supplied in the form of Ehrenfest- 
like relations. Electroproduction is modeled as a dis¬ 
sipative interaction of an energetic incoming electron 
with a medium. Since velocities of the incident particles 
are ultra-relativistic, a dissipative model for the Dirac 
fermion is required. The starting point of ODM are the 
Ehrenfest theorems for the coordinate ( x M ) and momen¬ 
tum (p^), which in the case of the free-particle Dirac 
equation (i.e., non-dissipative relativistic dynamics) read 

[22H2I1 

^V)=c< 7 v>, -jV)=o. a) 

ds ds 

Here, c is the speed of light, with p = 0,1,2, 3 are 
the Dirac matrices, and s is the proper time. In non¬ 
relativity the momentum and velocity operators are given 
by p k and p k /m , k = 1, 2, 3, respectively. However, be¬ 
cause a relativistic particle cannot move faster than light, 
the momentum and velocity operators have fundamen¬ 
tally different forms: p k and C 7 ° 7 fe , respectively. The 
former is unbounded, while the latter is bound by the 
speed of light [2Uj. Note that Eq. 0 is covariant. 
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The magnitude of a dissipative force of friction is 
proportional to the particle’s velocity relative to the 
medium, slowing the particle down. Therefore, the in¬ 
clusion of the friction leads to the following Ehrenfest 
theorems 

j t {x k ) = c( 7 V>, j t (p k ) = -2amc(7°7 fc ), (2) 


We require that these identities should be valid for any 
state p, hence 


op op 

~,dB sst .. , 

B'~— - t^-B = —Arnica . 

ox ox 


(7) 


where a > 0 is a friction coefficient phenomenologically 
characterizing the medium. Since the velocity is mea¬ 
sured with respect to the medium, Eqs. ([2| are non¬ 
covariant, and t is the time in the frame of reference 
attached to the medium. 

The Ehrenfest theorems |2]) can be reached from many 
other perspectives. For example, in an analogy with 
Ohm’s law, microscopically describing collisions [2B|, an 
Ohmic environment is coupled to the particle’s mo¬ 
mentum through the current, whose average for Dirac 
fermions is (7°7 fc ) 2? : , Mj. Additionally, a particle is 
coupled to the external photon field through the velocity 
operator, as prescribed by the minimal coupling in quan¬ 
tum electrodynamics [25] [25]. Moreover, a classical coun¬ 
terpart of Eq. ([2]), where all the operators are replaced 
by corresponding classical quantities, is closely related 
to the Langevin equation for the relativistic Ornstein- 
Uhlenbeck process mm, describing the velocity of a 
Brownian particle under the influence of friction. 

Equations ([ 2 ]) are consistent with a picture of open 
system dynamics. Hence, the average of an observable O 
is defined by 

(O) = Tr [Op] , (3) 

where p is a density matrix representing the state of the 
system, whose evolution equations is of the Lindbladian 
form 02HM] 

H = ca k pk + me 2 7 0 , 

D\p] = \ ( BpB t - \{p&B + , (4) 

where a k = 7 ° 7 fe and B is an unknown operator to be 
found following Ref. [2Tj . 

Hereinafter, we will consider one spatial dimension 
without loss of generality. Using the relations 

[B(x,p),p] = ihdB/dx , [B(x,p),x\ = —ihdB/dp, (5) 


that follow from the canonical commutation relation 
[x,p\ = ih and Eqs. and 0 , we obtain 


Tr 

Tr 
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dp 

dB 
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= 0, 


= —4imcTr [a 1 ^]. 


( 6 ) 


These are equations for the unknown operator B = 
B(x,p) defined as a function of the non-commutative 
variables within the Weyl calculus j55H57] , Perform¬ 
ing the Weyl transform, the operator equations Q turn 
into a system of equations for the complex-valued matrix 
function B(x,p) defined on the phase space 


B\x,p)* 

B\x,p)* 


dB(x,p) 

dp 

9B(x,p) 

dx 


dB^(x,p) 

dp 

dB^(x,p) 

dx 


* B(x,p) = 0, 

*B(x,p) = —Aimca 1 , 

( 8 ) 


where * denotes the Moyal product [551457] 

ih (^d~$ { d d^\ 

6XP 2 y dx dp dp dx J 

The system ([8]) has the following exact solution 

B{x) = —iVmcX — 2^/mcJX a x x, 


(9) 


( 10 ) 


where A is a constant with the dimension of length, whose 
physical meaning will be elucidated below. 

Even though D[p\ is non-translationally invariant, it 
can be replaced by the translationally invariant modifi¬ 
cation within the interval [—L,L\, 

D [p](x) ^ j dx 0 D[p\(x - x 0 ), (11) 

without changing the underlying Ehrenfest theorems ([2| . 

Let us find the non-relativistic limit of B(x). Note 
that B{x) is a matrix valued function of x. Calculating 
the eigenvalues of B(x), denoted by B e , and taking into 
account the shell-mass condition me = \JE 2 /c 2 — p 2 , we 
arrive at 


lim (B e e 5l7r / 4 ) = (iX + 2x signp)\/|p|/A, (12) 


Equation (121 coincides with the Lindbladian operator 


for the non-relativistic quantum friction recently ob¬ 
tained in Ref. [21], which leads to the interpretation of 
A as the thermal de Broglie wavelength, phenomenologi¬ 
cally characterizing the medium. 

Positron fraction in cosmic rays. Using a modifica¬ 
tion of the numerical method developed in Ref. [38], 
we propagate Eqs. 0. ( [To] ) , and ( |TT| ) in the phase 


space, where the density matrix p is represented by 
the Wigner function W(t,x,p). To construct an initial 
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state W{ 0, x,p) composed only of electrons, we apply the 
phase-space generalization of the projection operator j.'IDI 
to remove positrons from a Dirac gaussian spinor with 
the mean momentum po and the spatial standard devi¬ 
ation a. Time evolution generates positrons. To sepa¬ 
rate the state of created positrons W p (t,x,p) from the 
total Wigner function W(t,x,p), we project the elec¬ 
trons out from W(t,x,p), as prescribed in Ref. [59] . 
The fraction of positrons at the final time T is defined 
as N f W p (T,x,p)dx/ f W(T,x,p)dx, with N being a 
normalization constant. In the ultra-relativistic regime, 
which is employed here, p equals to the kinetic energy. 
For the propagator not to run into numerical stability 
issues due to a wide range of energies spanning the hor¬ 
izontal axis in Fig. |T| we scale the momentum p —> ep. 
The values of A = 400 and T = 7.68, in natural units, are 
chosen for numerical convenience and not varied during 
the fitting procedure. 

The five parameters a = 1 /a/3T4, p 0 = 4.8, er = 0.37, 
e = 22.5891, and N = 0.154 were used to fit the exper¬ 
imental points in Fig. |T| It is important to note that 
the overall shape of the curve in Fig. |T] is maintained 
for a wide range of values of the parameters. The first 
two (a and po) define the distribution of the incident 
electrons, the friction coefficient <r characterizes the in¬ 
terstellar medium, and the energy scaling constant e is 
introduced for a technical purpose described above. The 
position of the positron fraction minimum in Fig. |T| de¬ 
pends on po and cr; whereas, the slope of the distribution 
at higher energies is controlled by a. 

The developed Lindbladian model of electroproduction 
predicts that the positron fraction reaches a plato at high 
energies. The plateau seen in Fig. |T] is attributed to the 
fact that the friction force, causing electroproduction, is 
bounded because its magnitude is proportional to the 
velocity limited by the speed of light. Thus, pair produc¬ 
tion saturates for sufficiently high energies. Note that 
the experimental data in Fig. [T] cannot be fitted if the 
Ehrenfest theorems |2| are replaced by 

f t (x k ) = c( 7 V), j t (p k ) = -2 a(p k ), (13) 


where the dissipative force is proportional to the momen¬ 
tum, which is unbounded. 


Conclusions. In this Article we develop the master 
equation for electroproduction, a higher-order quantum 
electrodynamic process resulting in creation of electron- 
positron pairs upon ultra-relativistic electrons interact¬ 
ing with the medium. This model, being in an excellent 
agreement with observed data j7UT0| (see Fig. [I]), ac¬ 
counts for the unexpected rise of the positron fraction 
in cosmic rays, and thus, renders explanations involving 
dark matter decay unnecessary. According to our view, 
additional positrons are generated from cosmic ray elec¬ 
trons colliding with the interstellar medium leading to 
electron-positron pair production. 


Moreover, the ideas presented in the Article serve as a 
major step towards computationally clean formulation of 
quantum electrodynamics, which is bound to open new 
horizons in many-body quantum electrodynamics calcu¬ 
lations, where a discrepancy between theory and exper¬ 
iment has been observed HB HU. It has been known 
for a long time [42] that quantum electrodynamical ef¬ 
fects can be viewed as emerging from a particle inter¬ 
acting with a noisy environment, modeling the vacuum. 
Such a point of view currently enjoys a revived interest 
due to the need for computationally accessible alterna¬ 
tives for incorporating quantum electrodynamical effects 
into chemical structure calculations S33] - All these ap¬ 
proaches, nevertheless, have a limited rage of applica¬ 
bility because they represent ad hoc modification of the 
Dirac equation, describing dynamics of closed systems. 
A Lindbladian master equation can give a universally 
applicable and inconsistency-free description of a quan¬ 
tum system coupled to an environment. In this regard, 
the presented derivation of the master equation for elec¬ 
troproduction lays out a roadmap to accomplish compu¬ 
tationally efficient Lindbladian description of quantum 
electrodynamics. 
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